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The electroweak phase transition in the Two-Higgs-Doublet Mo- 
del is investigated. The Gibbs potential at finite temperature is 
computed with regard for the one-loop plus ring diagram contribu- 
tions. The strong first-order phase transition satisfying Sakharov's 
baryogenesis conditions is determined for the values of scalar field 
masses allowed by experimental data. The relation between the 
model parameters supplying the phase transition to be of the first 
order is derived. It is shown that a sequence of phase transitions 
is also possible. A comparison with results of other authors is done. 



1. Introduction 

Nowadays it is well known that the Minimal Standard 
Model (MSM) of elementary particles cannot explain the 
baryon asymmetry observed in the Universe. The mech- 
anism of generation of this asymmetry from the initially 
symmetric state was proposed by A.D. Sakharov ([l||, 
see also review and today is formulated as three 
baryogenesis conditions: 

1. Baryon number non-conservation. 

2. C- and CP-symmetry violation. 

3. Deviation from thermal equilibrium. 

In electroweak theory a deviation from the thermal 
equilibrium can be provided by the electroweak phase 
transition (EPT). The investigations showed that 

the EPT in the MSM is strong enough for the Higgs 
boson mass values that are incompatible with the mod- 
ern experimental bound m/j > 114.4 GeV. Monte Carlo 
simulations were used to study the EPT in the effec- 
tive three-dimensional gauge theory [4]. It was shown 
that a critical point exists in this model. For the expe- 
rimentally allowed rrih values the EPT becomes of the 
second order. In (Hj the EPT in the MSM in presence 
of external magnetic and hypermagnetic fields was in- 
vestigated. It was concluded that the third baryogenesis 
condition is not fulfilled, and the EPT becomes of the 
second order in strong fields. Thus, a deviation from the 



thermal equilibrium in the MSM is not strong enough, 
and Sakharov's baryogenesis scenario is not realized. 

In connection with this, the investigation of the EPT 
in extensions of the Standard Model is of substantial in- 
terest. One of the extensions is the Two-Higgs-Doublet 
Model (THDM) 0,0|. The THDM predicts four addi- 
tional scalar particles: a neutral particle iJ, a pair of 
charged fields H^, and a pseudoscalar particle Aq. As 
compared with the MSM, the THDM contains more free 
parameters in the scalar sector. 

The EPT in the THDM was investigated in Q. The 
obtained results were revised in fol by using an improved 
approximation (avoidance of the high-temperature ex- 
pansion, inclusion of the ring diagrams, and taking the 
experimentally measured mass of the top quark into ac- 
count). It was concluded that the third baryogenesis 
condition is fulfilled in the THDM for some parame- 
ter values. In [10.] a connection between the strength 
of the EPT and the one-loop correction to the triple 
self-coupling of the lightest Higgs boson was considered, 
and the deviation of the coupling value from the SM 
prediction was estimated. 

In [il|,[i3, as well as in the recent paper [isj, possible 
scenarios of the model behavior during the cooling were 
investigated. It was shown that the realization of a cer- 
tain scenario strongly depends on the parameter values 
of the scalar potential. The authors have computed the 
Gibbs potential in the THDM in a simple approximation 
considering only the first non-trivial finite-temperature 
corrections to the tree-level potential. Within this po- 
tential they observed a few interesting phenomena: 

1. A sequence of phase transitions (a second-order 
EPT breaks the electroweak symmetry; then a first-order 
EPT occurs). 

2. Exotic CP-breaking or charge-breaking minima can 
be realized during the cooling. 

The main goal of the present paper is to determine the 
domain of the THDM parameters, for which the strong 
first-order phase transition happens and the third baryo- 
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genesis condition holds. We use the approximation of 
thermal equilibrium. Following we compute the 

Gibbs potential in the one-loop order and include con- 
tributions of the ring diagrams. This consistent approx- 
imation, ensuring the minima of the effective potential 
to be real, allows us to check the results of fl^. 

We answer the question which is the domain in the 
space of model parameters that ensures the strong first- 
order phase transition. To estimate that, we introduce a 
certain relation between the model parameters and show 
that if this relation is satisfied, the system undergoes a 
strong first-order phase transition. 

In ^8i-,lQi] the special restrictions on the THDM pa- 
rameters were imposed, and the authors assumed that 
these restrictions ensure the symmetry breaking along 
the tan/3 = vijvx = 1 direction. Here, Vi are the vac- 
uum expectation values of the doublets. This was done 
to avoid difficulties with the two-stage phase transition 
lll . ll2i] . We do not restrict our investigation to this spe- 
cific case. Below we show that the noted assumption is 
not true in general, and a sequence of phase transitions 
may occur. The phase transitions happening may be ei- 
ther of the first- or second-order. In any case, the jump 
of the order parameter may differ essentially from that 
observed in [9]. 

The paper is organized as follows. Section [5] contains 
the necessary information on the Lagrangian and the pa- 
rametrization we use. Section |3] is devoted to the com- 
putation of the Gibbs potential. In Section |4] we present 
the obtained results on the phase transition. We also 
discuss the results obtained in |,S.il3i] and present the re- 
lation between parameters ensuring the first-order phase 
transition. Concluding remarks are given in Section [SJ 
The appendices contain a special information used in the 
main text. 



symmetry [rlfl^. 
1 



(2) 



This means the invariance with respect to the transfor- 
mation 



-ifl, IP2 ^ ^2- 



The neutral scalar fields Ci acquire the non-zero vac- 
uum expectation values (VEV) Vi and break the SU{2) x 
Uy{^) symmetry giving masses to gauge bosons and 
fermions. The mass mixing of the a^, Ci, and di fields 
takes place. The spectrum of physical particles is obtai- 
ned by the substitution 

~ af cos 7 + sin 7, = —a'^ sin 7 -f- cos 7, 
h — ci cos a + C2 sin a, H = —ci sin a + C2 cos a, 
X3 — di cos S + d2 sin S, Aq — — di sin S + d2 cos S, (3) 

where a-nd ^3 are the Goldstone modes. The expres- 
sions for a, 7, and 6 and their dependence on 111.2 are 
adduced in Appendix A2. 

In what follows, vi^2 denote arbitrary shifts of ci^2 
fields. These shifts define a "jump" of the order parame- 
ter during the phase transition. The extremum points of 
the potential are denoted as uoi,2- They are defined as 



dV 

dc. 



0. 



(4) 



2. Lagrangian 

The THDM Lagrangian differs from the MSM one in the 
scalar and Yukawa sectors. It can be written as 



The scalar sector is 
2 



Yuk 



;aue;c fixine+Khost ■ 



1 



E 



^V, (1) 



where y is a scalar potential. To simplify the analy- 
sis we restrict our consideration to CP-conserving vacua 
only. We consider the potential which possesses the Z2 



We consider the case where the VEVs are non-zero for 
both doublets, i.e. tioi,2 7^ (see Appendix Al). 
The gauge-boson masses are 



2 ff^ 2 
1^^ = . 



In di ■ [l3| a similar scalar sector is considered. The 
main difference is that there are the additional terms 
fi^(p\(p2+ h.c. softly violating the Z2 symmetry. As was 
noted in [8|, the influence of these terms on the EPT 
strength is small (though their presence allows one to 
introduce an additional CP violation). Another distinc- 
tion is the parametrization of the scalar field couplings 
(see Appendix A3). 
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As we mentioned above, the analysis in (^,[^ was sim- 
plified by restricting the possible tree-level parameter 
values to 



2 _ 2 \ _ \ 



(5) 



In what follows we consider this case separately. 

The general parametrization for the Yukawa interac- 
tion is 



^Yuk — — 



z 



/l *=i 



+ Gu, [flmu)n + ifu)n^ffi] } , 

f _ 1-75 f fu\ r _ 1+75 f 



(6) 



Here, ip^ = i<J2(pl- The heavy quarks - t and b - only 
are of importance for the Gibbs potential, so we neglect 
the CKM mixing. 

This Yukawa Lagrangian leads to the existence of 
flavor-changing neutral currents (FCNC) |l7|. According 
to the Glashow- Weinberg theorem [l^, the dangerous 
processes with the FCNCs can be excluded at the tree- 
level if all fermions of a given electric charge couple to 
no more than one Higgs doublet. The most popular pa- 
rametrizations which respect these restrictions are the 
THDM type I and the THDM type II 



• in the THDM type I all fermions are decoupled 
from the second doublet, i.e. Gd,2 = Gu,2 = 0; 

• in the THDM type II the u, c, and t quarks couple 
to the first doublet, while d, s, and b couple to the 
second doublet, i.e. Gd.i — Gu.2 = 0. 

We consider the THDM type II parametrization. The 
main reason is that it represents a low-energy limit of 
the Minimal Supersymmetric Standard Model. 

There is an interesting feature in the THDM type II. 
It follows from the expressions for the quark masses at 
the tree-level minimum of the potential: 



mtfi^Gt,iVoi, nibfi ^ Gb,2Vo2- 



(7) 



Let vq2 be one or two orders of magnitude smaller than 
vqi. Then Gt,i and Gb,2 couplings have to be of the 
same order of magnitude to preserve the quark mass ra- 
tio. This is an essential difference from the MSM and 
the THDM type I. By noting this, we find that in the 
THDM type II case it is necessary to include the b quark 
contributions to the Gibbs potential and Debye masses. 



To compare our results with that of [9] we also consider 
the case of Gt_i — Gt,2- In this parametrization the 
contribution of the b quark to the Gibbs potential is 
also negligibly small. We will refer to this together with 
restrictions ([5]) as the doublet-universal parametrization. 

In the present paper, all calculations are carried out in 
the Feynman-'t Hooft gauge. The gauge-fixing functions 
are 



G 







(8) 



Then, the gauge-fixing part of the Lagrangian reads 

^gaugc fixing ^ ^ 9 I X^ G + G j . (9) 
^ \a=l J 



The quadratic terms in the Faddeev-Popov sector are 



host 



(10) 



where ^ is a gauge-fixing parameter. For arbitrary ^, the 
gauge-boson propagator is 



iD^^'ip) = 5 



9^" + - 1) 



The fermion and gauge field sectors of the Lagrangian 
are taken to be 

/l 

/r 

C.^^^F'^'^F.^-^FrF^^. (11) 
We proceed with the calculation of the Gibbs potential. 



3. Gibbs Potential 

The calculation of the Gibbs potential was discussed 
in numerous papers [13, Ell (see also review [l^). We 
compute the effective potential in the following standard 
way: 



(12) 
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where Vtrcc is the tree-level potential, V^^^ is the one-loop 
correction at zero temperature, V^i is the one-loop finite- 
temperature correction. The tree-level part is obtained 
by substituting ipQi into ([2]) and reads 



Vireo(Wi) = --^{PlVi + ^2^2) + ''^l^l + ^2^2 + 
+ {X3 + X5Hvl 

Now, we consider V^^^ and V^. 



(13) 



3.1. One-loop contributions at zero temperature 

The regularized contribution of a field of mass m to Vg 

is E3 



-fm4(ln(W)-3/2 + 7-^; 



(14) 



We use Schwinger's proper time regularization with the 
regularization parameter sq. It has to be set to zero at 
the end of calculations. 

The general expression for the scalar field mass is 



±JiCivf + C2VI + C3)2 + iDiViV2) 



(15) 



Here, Bi, Ci, and Di are some combinations of the tree- 
level VEVs and couplings. There are four pairs of the 
scalar fields: h and H, and H^, xs and Aq. The 
sign in Eq. (|15p . "— " or "-|-", corresponds to the mass of 
one field of a pair. For example, we have the "— " sign in 
case of the X3 mass and the sign for the Aq mass (see 
Appendix A2). In the sum of both field contributions, 
the term 



1 



- / 4 



V m_ 



(16) 



appears. It is cancelled out by the term coming 
from a high-temperature expansion, when the finite- 
temperature corrections are taken into account. How- 
ever, we do not use this expansion in our calculations, 
and, therefore, the explicit cancellation does not occur. 
This term results in cumbersome quantum corrections 
to Vg. 

The contribution coming from fermions, gauge bosons 
and ghosts is given by (fT4|) with regard for the factor A, 



A = 



-1x4x3, 
-2, 

3+e^ 



quark, 
ghost, 

gauge boson. 



(17) 



This factor accounts for the number of degrees of free- 
dom and the color states of fields. 

We choose the renormalization conditions preserving 
the tree-level vacuum energy value, VEVs, and mass 
terms. They are taken to be 



K(wOi) = VtrccivOi), 



0, 



dvf 



+ 12Aiwgi + 2(A3 + A5)wo 



dvidv2 



0. 



12A2W^2 + 2(A3 + A5)w^i, 



(18) 



Since the renormalized contributions from the scalar 
sector are cumbersome, we do not adduce them here. 
They could be obtained easily by using a symbolic cal- 
culation software. The renormalized contributions of a 
fermion, a gauge boson, or a ghost field read 



vy{m) = 



A 
647r2 

2 \2 



'(ln(- 



(19) 



where mvac is the field mass value at Vi = voi. 

The complete temperature-independent part of the 
Gibbs potential is 



-1-2(3 + e)V,'^'{mw) + (3 + eW^^'^imzy 



(20) 



where V}'!: „, V^'^+ „+ and V}'^ a are the scalar field 
contributions. 

3.2. One-loop contributions at finite 
temperatures 

Finite-temperature corrections are calculated by using 
the Matsubara formalism. For the contribution of one 
bosonic degree of freedom we have [Tsj 



4 ^ 

V^{m) J dxx^ In (l - cxp{-^/ + P^m^)'^ , 



(21) 
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where /3 is the inverse temperature. We recah that Mat- 
subara's frequencies for ghost fields are even [20J. For 
one fermionic degree of freedom, we have 



CO 

T4(m) = — / dxx^ In \\ + exp(- Va^^ + Z^^m^) 



(22) 



The degree-of-freedom factors for fermions and ghosts 
stand in ([T7)) . The contribution of the massive gauge 
field is 

ys'^"s<=(m) = 3T/^(m) + (23) 

The last term in cancels a part of the ghost field 
contribution. 

3.3. Ring diagram contributions 

As is well known, an imaginary part of the one-loop 
Gibbs potential arises at small Vi. It comes from the 
scalar sector contributions at finite temperature and in- 
dicates the instability of the system. Gauge bosons are 
massless in the symmetric phase that leads to infrared 
divergences. These shortcomings of the one-loop effec- 
tive potential can be avoide d by adding the ring-diagram 
(Fig. [T]) contributions [illilHiil. These dia grams in- 
troduce additional finite-temperature corrections to the 
masses of bosons that result in the terms of the order 
^ g'^ or ~ Aj m Vg- 

The ring-improved finite-temperature correction to 
the Gibbs potential in the scalar field case is 

V^(m(r))==^ Jdxx^ In (l~exp{-^/x^+l3^m^{T))j , 



m2(r) = +Sm'^(T), 



(24) 



where 5m? {T) denotes the Debye mass of a field. This 
correction ensures that the imaginary part for the finite- 
temperature contribution is absent. The Debye mass is 
defined through the polarization tensor 11 of a field taken 
in the infrared limit (23| 



5m^{T) = n(/co = 0,fc 0), 



(25) 



where k is the four-momentum of a field. 

For gauge fields, the Debye mass is defined as — IIoq 
in the infrared limit. The ring-diagram contribution to 
Vq from each massive gauge boson is Q 



(26) 



Fig. 1. Ring-diagram contributions for scalar field. Black blobs 
denote Debye masses 

The one-loop Debye masses of the Higgs fields in the 
THDM type II are 

5ml{T) = U2Ai + ^G^) cos^ a + (2A2 + 



+ -Gl) sin^ a 



A3 + A4 + 4A5 3.g2 + g- 



12 



6 



16 



Smjj{T) = (^(2Ai + -G^) sin^ a + (2A2+ 

^1^2^ 2 , A3 + A4+4A5 , 3g^+g'^ 
+ -G,) cos a + + 



(27) 



For other scalar fields, the Debye masses are given by 
similar expressions. The difference is that the angle a is 
replaced by 7, (5 from Eq. ([3]). 

The complete Higgs-sector contribution to the finite- 
temperature part of Vg is 

V^iv,) ^ V^imhiT)) + V^imniT)) + 2V^im^± (r))+ 
+2VS{mH±{T)) + V^{m^,{T)) + V^{mA,{T)). (28) 

The gauge-boson Debye masses are 

6ml.{T)^2g-T\ ^^1^^) =^^^±^T- . (29) 

For Faddeev-Popov ghosts, the Debye mass is zero in 
the leading order in T. 

In Appendix B we give the Debye masses of all fields. 
We note that these corrections are ^-independent. 

The finite-temperature gauge-field and ghost contri- 
bution to Vq is 

V^{v{) = W^imw) + iV^{mz) - 2V^{^^mw)- 
-V^iVl^z) + 2F™s(^^) + l/"ng(™z). (30) 

For the final expression of Vq in the THDM we have 

Vaivi) = V^{v^)+V^{v,)+V^{v^)~l2V^{mt)-l2V^{mi,). 

(31) 

The minimum value of Vg is gauge-independent. The 
^-dependence is cancelled out between the gauge boson, 
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0.2 0.4 0.6 0.8 

vJT 

Fig. 2. Two distinct minima signaling the EPT of the first order 
are realized. The critical temperature is Tc = 125.65 GeV 



Goldstone field and ghost contributions. Physical quan- 
tities are also gauge- invariant (24i |. 

However, the Gibbs potential is gauge-dependent [l7| 
at arbitrary values of Vi. For numerical calculations one 
has to choose the value for ^. We set ^ = 1. In this 
case the Goldstone field masses equal to the masses of 
corresponding gauge bosons. 

4. Phase Transition 

The third Sakharov condition is fulfilled if the first-order 
EPT is realized, and the order parameter jump is greater 
than 1 



where Tc is the critical temperature. 

In the THDM there are two order parameters - vi/T 
and vijT . In general, several jumps of the order para- 
meters may occur, i.e. single phase transitions, as well 
as sequences of phase transitions, may happen. During 
a series of EPTs the system goes to an intermediate vac- 
uum state, where the symmetry is broken for one doublet 
only. We discuss possible scenarios of phase transitions 
in the THDM and find the domain in the parameter 
space, for which 5v/Tc is large. 

We consider Ai.2,3,4,5 and one of the VEVs voi,2 as the 
free parameters of the model (see Appendix Al). 



Possible scenarios 

The plots of 100(Vg(wi,W2) - Vg{Q,Q))/T'^ versus wi/T, 
V2/T are shown in the figures. The blue areas repre- 
sent lower values of the Gibbs potential. The parameter 
values for the figures can be found in Appendix C. 
First, we consider the THDM type II. 

1. Fig. [21 The form of Vq indicates a strong first-order 
EPT. The order parameter jump is 

5v = V^^pM = 1.01. 

For this set of parameter values, the tree-level scalar field 

masses are 

ruh = 119 GeV, tuh = 131 GeV, mjj± = 181 GeV, 
TOAo = 338 GeV. 

This scenario is the most favorable for successful baryo- 
genesis. 

2. Fig. [3l The sequence of phase transitions is gener- 
ated. The weak first-order EPT breaking the symmetry 
along the directions tan /3 = or tan /3 = -\-oo happens 
(the former case is shown in Fig. 3, a). Then next weak 
first-order EPT follows (Fig. 3,&). The Gibbs potential 
minimum is now located along the < tan/3 < -|-oo di- 
rection in the (wi, V2) plane. The scalar field masses in 
the tree-level approximation are 

ruh = 114 GeV, mn = 132 GeV, to^± = 181 GeV, 
TOAo = 266 GeV. 

The successful baryogenesis cannot be realized. The 
thermal equilibrium approximation used could be unre- 
liable. This is because of non-equilibrium processes hap- 
pening after the first phase transition. In fact, this series 
of phase transitions, as concerns its consequences, could 
substitute one strong enough first-order phase transition. 
The calculation of characteristics for a phase transition 
of such type requires other methods and additional in- 
vestigations. 

3. Fig. [31 The sequence of phase transitions happens. 
A strong first-order EPT breaking of the symmetry along 
the tan /3 = or tan/? = +00 direction is generated (the 
former case is shown in the figure) . Then the system un- 
dergoes a second-order phase transition, the tan (3 value 
being finite and non-zero. The scalar field masses are 

TO,, = 120 GeV, mu = 201 GeV, to^± = 322 GeV, 
TOAo = 429 GeV. 

This scenario is acceptable for the successful baryo- 
genesis. 
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0.05 



0.02 




0.012 



0.008 



0.002 




0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.6 
b 

Fig. 3. Sequence of the first-order phase transitions. The critical temperatures are (a) Tc = 128.18 GeV, (b) Tc = 127.66 GeV 



4. A sequence of second-order phase transitions is pos- 
sible. Such scenario is reahzed, for example, if the scalar 
field masses at the tree level are 

m,, = 114 GeV, mn = 162 GeV, m^± = 181 GeV, 
ruAo = 276 GeV. 

In this case the baryogenesis is not realized. 

These results also hold for the THDM type I qualita- 
tively and, in most cases, quantitatively. The main dif- 
ference from the THDM type II is that the b quark influ- 
ence on the Gibbs potential is much weaker. Therefore, 
a quantitatively different picture should be expected for 
small tan/3o = vo2/vai values, namely tan/3o < 0.1. For 
this domain of the parameter space the b quark contri- 
bution to Vg in the THDM type II is non-negligible. 

Note that the parameter values considered correspond 
to Sector I of the parameter space defined in fl3|. In this 
sector, the authors observed a second-order EPT. They 
proposed that a more detailed investigation based on 
a consistent corrected effective potential could predict 
the EPT of the first order. This is just what we have 
observed. 

4-2. The doublet-universal parametrization 

Since there is a possibility of a sequence of phase tran- 
sitions, it is interesting to discuss the assumption about 
the symmetry breaking made in [8l-[lo|. In [oj, a strong 
first-order EPT was observed. It was assumed that the 



0.12 




- 1.2 



0.04 



vJT 

Fig. 4. The first-order EPT in the sequence of phase transitions. 
The critical temperature is Tc = 153.3 GeV 

symmetry breaking happens along the tan /? = 1 direc- 
tion, if restrictions ([S]) are imposed. This allows one to 
express all free parameters of the potential in terms of 
the tree-level masses of fields. The domain of the mo- 
del parameter space considered in corresponds to the 
condition for scalar field masses 



(33) 



In [9] the tree-level masses of scalar fields were used as 
free parameters instead of the couplings. However, the 
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mass values do not define uniquely the values of cou- 
plings. Really, let us consider the h field mass case. The 
tree-level expression of to/j at the minimum point of the 
potential is 

ml = 4AiUoi -I-4A2W02- 
-4 



By applying ([5|), we obtain 
ml=Avl^ (2Ai- IA3 + A5I). 



(34) 



Since A3 < (see Appendix A2), there exist two pos- 
sible values of A5 corresponding to the same value of 
m/i, namely 



A± = -A3 ± (2Ai - 



mr 



■)■ 



(35) 



It appears that, at A5 — Min(A5 , X'^) for the parameter 
values resulting in (1331) . the symmetry is broken along 
the tan /3 = 1 direction, as it was assumed. However, if 
A5 is taken to be Max{X^ , A^), then the evolution of the 
system is completely different. The system undergoes 
the sequence of first-order phase transitions. The min- 
ima with broken symmetry appear along the tan (3 — 
and tan/3 ~ +00 directions. After that, the Gibbs po- 
tential develops another minimum along the tan/3 = 1 
direction. This sequence is shown in Fig. [5l We can 
conclude that the assumption mentioned is not always 



true. These model parameters correspond to the scalar 
field mass values 

ruh = 120 GeV, ttih = mH± = ruAf, ~ 250 GeV. 

Note also that in this parametrization three degener- 
ate vacuum states may coexist. 

4-3. Relation between model parameters 

It can be seen from Appendix C that a small change of 
any parameter value may result in a significant change in 
the system's evolution (compare the scenarios shown in 
Figs. [5] and [3]). So it is interesting to determine any re- 
lation between model parameters which provides a large 
Sv/Tc value, and, hence, a strong first-order phase tran- 
sition. Let us turn to this problem. 

From geometric reasons it is natural to assume that a 
jump of the order parameter is large enough if the non- 
trivial vacuum appears in the < tan /3 < -l-oo direc- 
tion. This happens if the symmetry is broken for both 
Higgs doublets simultaneously. The analytical form of 
this condition reads 



M?(r) = o, 
Mi(r) = o. 



(36) 



The values of couplings ensuring this condition can be 
found by using the high-temperature expansion of (|21l) 
and (|22l). 

Since the quantum corrections are cumbersome in the 
THDM, we discuss the application of ([36|) considering a 
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toy model as an example. Then we apply this relation 
to the THDM type II. 

The toy model has to possess two main properties of 
the THDM - the spontaneous symmetry breaking and 
the mass mixing of scalar fields. The Lagrangian we use 
is 



C 



V, 



1 



y = -^(pI^I + f^l^l) + ^i^i + ^2^t + h^lvl (37) 

where (pi are real scalar fields. The Lagrangian is invari- 
ant with respect to the transformation (pi — > —'Pi- 

Let us shift ipi by arbitrary values Vi and obtain the 
mass eigenstates as those in Eq. ([3]). We denote a new 
pair of fields as /ii_2- The tree-level VEVs wpi are obtai- 
ned from the relations 

-^ll+A\lvl + 2\^vl = 0, 
-^l+AMvl + 2A3W^ = 0. 

The tree-level potential is 

T4,ec = + M2«2) + + X2V2 + A3W?v|. (38) 

The one-loop corrections to the Gibbs potential are given 
by Eqs. (fH|) and (^1]) . We use the MS renormalization 
scheme. Then the one-loop vacuum corrections can be 
written as 

The high-temperature expansion of (|211) looks as follows: 

, tt'^T^ T^m^ Trrv' 
Vrp im) 

1 



90 



4n 



24 
5.41). 



127r 



For the one-loop Gibbs potential, we have 



+V^{mh,) + V^{mn,). (39) 
The functions ^i.2{T) are the factors at 2 in (f39| : 

^xl{T) = ^4 + -^{l2X2fil + 2X3fil)x 

X hn + 10.82j - ^(12Ai + 2A3), 
1 



M^(r) = ^4 + ^(i2AiA.^ + 2A3Aif); 

rpi \ rp2 

12 




I 



0.4 



Fig. 6. The phase transition in the THDM type II with the optimal 
parameter values. The critical temperature is Tc = 125.39 GeV 

There are seven parameters in the toy model - Ai.2,3, 
/ii 2, Wo 1,2- Four of them are free. Note that we took the 
condition v^i -\- Vq2 = const into account. The way of 



X ( In — + 10.82 ) - — (I2A2 + 2A3) 



(40) 



using (p6| is the following. We take arbitrary numerical 
values for any three parameters (for example, Ai^2.3), and 
the value of the remaining parameter (wqi) is determined 
by solving the system (p6|) . 

The fulfillment of (|36p is sufficient for the symmetry 
to be broken along the < tan/3 < -l-oo direction. 

4.4. THDM case 

Let us apply the procedure described above to the 
THDM type II. In there are nine parameters - 
Ai,2,3,4,5, Ml. 2, ^^01, 2- Six of them are free (we choose 
Ai,2,3,4,5 and Vol). We set the values of the couplings 
(listed in Appendix C), then determine ^i^2{T), and 
solve ([36|) to obtain the uoi value for each set of the pa- 
rameters. We will refer to the parameters obtained this 
way as the optimal parameters. The Gibbs potential 
with the optimal parameters taken is shown in Fig. [6] 
The Ai^2,3,4,5 values are the same as those in the case 
shown in Fig. [21 and the jump of the order parameter is 
somewhat larger: 5v/Tc = 1.02. 

If the parameters are optimal, and the fields Aq, 
are heavy enough {mu± > 80 GeV, [25J), the strong 
first-order EFT is realized in the THDM. For smah devi- 
ations from the optimal parameters, the third Sakharov 
condition is still fulfilled. As a result, we see that the 
baryogenesis condition is satisfied in a wide domain of 
the parameter space. Five of the six free parameters of 
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the model can be set to the values which are constrained 
by modern experimental bounds and stability require- 
ments for the tree-level potential. Then the remaining 
parameter can be calculated by using ([36]). 

5. Conclusions 

In the present paper we have investigated the electro- 
weak phase transition in the Two-Higgs Doublet Model 
on the base of the ring-improved one-loop Gibbs poten- 
tial. We found that there is a wide domain in the para- 
meter space of the model, for which the third Sakharov's 
baryogenesis condition is fulfilled. The values of the pa- 
rameters entering the potential correspond to the scalar 
field masses that are compatible with modern experi- 
mental constraints. The parameter values for this do- 
main can be found by using the introduced relation p6p . 

The EPT kind depends strongly on the tree-level pa- 
rameter values. We have observed that the single phase 
transitions, as well as the sequences of transitions of the 
first and second orders, may happen. 

We have concluded from our analysis that the restric- 
tions fif — fj.2j Xi — A2 do not ensure that the symmetry 
breaking is realized along the tan/3 — 1 direction, as it 
was proposed in . We have seen also that in the mo- 
del studied in [8,9], three degenerate vacua may coexist. 
Some of them can be realized as overcooled states. 

Nowadays, there are few essential experimental con- 
straints on the THDM parameter values: the lower 
bound on the mass of a charged scalar field in the MSSM 
and the bounds on tan /3 (see review in [25|). In this 
situation the parameter values, for which the successful 
baryogenesis is possible, could be used as certain refer- 
ence points in the study of a model extending the MSM. 

APPENDIX A 

Information on the scalar sector 

In this appendix we present the constraints on parameters of the 
scalar sector, the expressions for scalar field masses and the rela- 
tions between parameters in different parametrizations. 

1. Tree-level potential properties 

Six of the nine scalar sector parameters are free due to the mini- 
mum conditions for the tree potential 

-fii + 4\iv^i + 2(A3 + A5)i>g2 = 0, 

-til + 4A2fg2 + 2(A3 + X5)v^i = 0, (Al) 
and the VEV vq is known: 

''01+^02 =^0 = (246GeV)2. 



For the stability of the tree-level vacuum, the potential value at 
large Vi values has to be positive. This translates into Sylvester's 
criterion for the quadratic form 

Xivf + A2t>| -I- (A3 + \5)vlvl. 

Then the scalar field couplings are restricted to the conditions 



Al > 0, A2 > 0, 4A1A2 > (A3 + As)^ 



(A2) 



In the present paper we consider the case where the potential 
minimum is realized at uqi 0, V02 7^ 0. This is because other 
cases Doi 7^ Oj ^02 = or vice versa can be reduced to the MSM 
case with several additional fields, and no mass mixing is present. 
As it was shown in |ll|jl2| , the latter case also leads to a possibility 
of successful baryogenesis through a two-stage phase transition. 

2. Scalar field masses 

From ||3}, we derive the expressions for the angles a, 7, and S. 
The a angle and the h, H masses are 

A3 



tan2Q! : 



A2 - Al 



"ift, H = Al -I- A2 T V (A2 - Ai)2 + Aj, 
Al = + 6Aii)2 + (A3 + X5)v1, 

A2 = -^^J.^+ 6A2i>i + (A3 + X5)v1, 
A3 = AviV2(\3 + A5). 

The S angle and the X3, Aq masses read 
B3 



(A3) 



tan 25 ■ 



B2-B1 



B1+B2T J(B2- Bi)2 +Bl, 



1 



Bi = --tii + (2X1 + i{g' + g'^)/8)vt + {Xi + \5)vi, 

B2 = -\y'l + (2A2 + £.(9' + a'^)/?.)vl + (A4 + Xr,)vl, 

B3 = 2t,it;2(A3 - A4 -f ^(g^ + 9'^)/8). (A4) 
The 7 angle and the x^i masses are calculated to be 
2C3 



tan 27 : 



C2 — Ci 



m' 



C1+C2T y(C2 -Ci)2 +4Ci 



X±,H± 

Ci = 'l^^l + {2Ai + £g^/8)vf + (A4 + X5)vl 

C2 = -\i4 + (2A2 + £.g''/?>)vl + (A4 + X5)vl 

C3=viV2{Xa+£g^/8). (A5) 
At Vi = VQi, the Aq and masses are gauge-invariant, and 
the Goldstone fields masses are \/^mz, v^™W- The masses of the 
h, H, Ao, and H± fields are 

ml = 4Ai»;gi +4A2t)o2- 

-'^^J{X2vf,2 - Ai»^gi)2 + (doi1'02(A3 + A5))2, 
mjj = 4Ai?;gi +4A2»^o2-l- 
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™Ao = 2(A4 - A3)(l)gi + VI2), 

ni]j±=-2\3{vl^+vl^). (A6) 

These expressions yield the following restriction on the scalar field 
couplings: 

A3 < 0, A4 - A3 > 0. (A7) 

These constraints ensure that the scalar fields Aq and are 
physical ones. 

3. Different parametrizations 

The expression for the THDM potential used in fis] is 
V = —- ^mliXi + m222-'2 + mfjC^Js + 2^3)] + 



+ ' + A32:i2:2 + A4X34 + ^ ^ % 

where 



(A8) 



t t t WO 

Xl=(^l(/Pl, X2=ip2V>2, Xs=ip[ip2, ipi = —=\ 

v2 V ^« 

and we consider the CP-conserving case only. By comparing this 
expression with (Ell and lfT3)l . we obtain the relations between the 
parameter values in different parametrizations: 

mf]^ = , = 2/i2i Ai = 8A1, A2 = 8A2, 

A3=4A5, A4 = 2(A3-A4), A5 = 2(A3 + A4). (A9) 

Note that the parameter m\2 is absent in the potential investigated 
in the present paper. 

APPENDIX B 
Debye masses 

Finite-temperature corrections to the h and H masses are given 
by 1 I27I I. The Debye masses of the x^i H^, xSi ^0 fields are 
calculated to be 

Sral^ = T2 ('{2Ai + ^G?) cos^ 7 + {2X2 + ^Gg) sin^ 7+ 



A3 + A4 + 4A5 3g2 -f^ g'- 



6 16 
Sm%± = (^{2X1 + ^G?) sin^ 7 + (2A2 + ^G^) cos^ ■y+ 
, A3 + A4 -I- 4A5 3g2 -I- 



6 



16 

1 
2 

A3 + A4 + 4A5 , 3g2+3'2 



<5m2^ = ( (2Ai + ^G?) cos^ S + (2X2 + ^Gg) sin^ S+ 



6 



16 



5m^^ = ( (2Ai ^G?) sin^ 5 + {2X2 + ^gI) cos^ 5+ 



+ 



Xz+X4+ 4A5 , -I- g>'- 



6 ■ 16 y ■ '"'^^'^ 

For the THDM type I and for the doublet-universal parametriza- 
tion, we have the following differences: 



1. THDM type I. The G5 coupling is always small as compared 
to the Gt coupling and can be neglected. 

2. The doublet-universal parametrization. In this case, Gf, can 
be neglected. Also, one has to make substitutions into IIBll l using 
the prescription 

Gj sin^ a, 7, (5 — ^ G^ (sin o, 7, 5 -|- cos a, 7, (5)^ , 
Gt cos^ a, 7, (5 — > Gt (sin o, 7, 5 — cos a, 7, 5)^ . 

We present the contributions to the gauge boson Debye mass 
from different sectors of the model. The Z boson Debye mass is 

^ 3 g^+g'^ ^ 12 32+^/2 ' 



2^2 



3 92 + g,2 

For the W boson mass, we have 
Sm%^ = 5m\^ + 5m{^ + 5m^ , 

The s, /, and g superscripts denote the contributions coming from 
the scalar, fermion, and gauge boson sectors, respectively. 

APPENDIX C 

Parameter values for figures 

In this appendix, we adduce the parameter values used in the 
calculations. The t and h mass values were taken 175 GeV and 4.2 
GeV, respectively. 

Table 1. The parameter values for used for 
plotting the Gibbs potential 



fig- 


Ai 


A2 


A3 


A4 


As 


f0l/fO2 


m 


0.045 


0.135 


-0.27 


0.675 


0.27 


1.91 


m 


0.045 


0.135 


-0.27 


0.315 


0.27 


2 


[3 


0.095 


0.2375 


-0.855 


0.655 


0.855 


2.64 


El 


0.16 


0.16 


-0.52 





0.72 


1 


1 


0.045 


0.135 


-0.27 


0.675 


0.27 


1.86 



Table 


2. The parameter 


values resulting in a sequence 


of second-order ph 


ase transitions 






Ai 


A2 


A3 


A4 


As 


l'0l/l'02 


0.045 


0.135 


-0.27 


0.675 


0.27 


1.22 
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